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Abstract: In some extra dimension theories with a TeV fundamental Planck scale, black 
holes could be produced in future collider experiments. Although cross sections can be 
large, measuring the model parameters is difficult due to the many theoretical uncertainties. 
Here we discuss those uncertainties and then we study the experimental characteristics of 
black hole production and decay at a typical detector using the ATLAS detector as a guide. 
We present a new technique for measuring the temperature of black holes that applies to 
many models. We apply this technique to a test case with four extra dimensions and, using 
an estimate of the parton- level production cross section error of 20%, determine the Planck 
mass to 15% and the number of extra dimensions to ±0.75. 
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1. Introduction 

In recent theories with extra dimensions the fundamental Planck mass, Mpl, can be as low 
as the TeV scale [jl|, B ||, Q making the trans-Planckian regime accessible for future high 
energy colliders. These models aroused great theoretical interest because they address the 
weak-Planck scale hierarchy. In these theories microscopic quantum black holes could be 
produced at energies higher than the Planck mass at the LHC || |6|, [7|, ||, ||, [H], |ll|, |l2|| . 
Once produced, the black hole would decay very rapidly to a spectrum of particles by 
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Hawking radiation. Assuming that all the Standard Model matter and gauge fields are 
confined to the physical three-branes in a higher dimensional space, it has been shown that 
most of the black hole decay products are Standard Model quanta emitted on the brane Q 
and are therefore visible experimentally as very spectacular events. 

We stress that quantum extra dimensional black holes in no way constitute any threat, 
being distinguished from the more familiar astrophysical variety by being much lighter and 
highly unstable. The astrophysical variety is much too heavy to be produced in current 
or planned collider experiments. Hereafter, all discussion of black holes relates only to the 
extra dimensional variety. 

In the large extra dimensions scenario (l|, |2| at distances small compared with the size 
of the extra dimensions and in the warped scenario ||, ||] at distances small compared to the 
curvature scale of the geometry associated with the extra dimensions, black holes with a 
horizon radius, rpHi smaller than the size of the extra dimensions can be treated as higher- 
dimensional objects located on the brane and extending along the extra dimensions. It has 
been shown that these small black holes have modified properties, e.g. they are larger and 
colder compared to a 4-dimensional black hole with exactly the same mass ||. 

The black hole discovery potential is critically dependent on the value of Mpl. Short 
scale gravity experiments and particle collider experiments provide limits on the funda- 
mental Planck scale. However for smaller values of the number of extra dimensions, n, 
the more stringent constraints come from astrophysical and cosmological data, albeit with 
larger uncertainties. It is widely agreed that the one large extra dimension scenario is ruled 
out by such data. The present collider limits 1 on the Planck scale range from 1.3 TeV for 
n = 2 to 0.3 TeV for n = 6 arising from the production of real (from missing energy signa- 
tures) or virtual Kaluza-Klein gravitons at the Tevatron Run I and LEP II |l3|, [l4|, |l5[ 16]. 



For a comprehensive recent review of these constraints see, for example, p7f[ . The LHC, 
with a centre of mass energy of 14 TeV, offers a good opportunity for black hole production 
if Mpl ~ TeV. The very large cross section for production of black holes not too much 
heavier than the fundamental Planck scale corresponds to a production rate of a few Hertz 
at the LHC design luminosity. 

In the following sections, the process of the black hole production and decay is reviewed 
(section ||), followed by a description of the CHARYBDIS |l8|, event generator (section |||) . 
We then present a review of the principal theoretical uncertainties (section before moving 
on to experimental discussions. The characteristics of black hole decays are presented in 
section [| followed by a discussion of the measurement of the black hole mass in section |6[ 
We then discuss ways of determining the Planck mass (section ^) and finally, in section [8|, 
we study methods of determining the number of extra dimensions. Throughout, we have 
used the ATLAS fast simulation software |2(J to give a description of a typical detector 
and we have used the full simulation [^] to verify the main results. 

2. Black hole production and decay 

In the black hole event generator CHARYBDIS, which has been used in these studies, the 



1 Limits are given in the convention of b) which is used throughout this paper (see section 
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black hole production is treated as a semi-classical process (black hole mass, Mbh S> Mpl) 
and it is assumed that the extra dimensions are large (3> tbh)- For black hole masses close 
to Mpl this semi-classical approximation is not valid and a theory of quantum gravity 
would be required to calculate the cross section. To be within the semi-classical domain 
we restrict the mass of the black hole to be Mbh > 5 Mpl- By geometrical arguments the 
semi-classical parton-level cross section for black hole production would be || 

a(s = M| H )«7rr BH (2.1) 

where V§ is the centre-of-mass energy of the colliding particles (see equation 2.4]). 
The radius for a Schwarzschild black hole is 

where we have used the convention Mp^ 2 = l/G( n+4 ) where G( n+4 ) denotes the n + 4 
dimensional Newton's constant ||], so for a fixed black hole mass, the cross section is lower 
for a higher number of extra dimensions. This convention has been used throughout this 
paper. 

The decay of a spinning black hole comprises three phases H: 1) the balding phase, 
in which the black hole loses its 'hair' (associated with the multipole moments) by the 
emission of radiation; 2) the evaporation phase, which starts with a brief spin-down phase, 
shedding away its angular momentum, followed by the Schwarzschild phase, emitting a 
large number of quanta which reduce the mass of the black hole; 3) finally the Planck 
phase (also called the remnant decay), when Mbh approaches Mpl, in which the final 
decay takes place by the emission of a few quanta. 

CHARYBDIS only models the Schwarzschild phase which is expected [||] to account for 
the greatest proportion of the mass loss. The energy spectrum of decay products is ap- 



proximately black body with corrective 'grey-body' factors [19], 7, which the generator 
includes. The spectrum for a fixed temperature black hole is 

dN E 2 -f 

dE ^ ( e £/TH T 1) T£+ 6 (2 ' 3) 

the denominator includes a spin-statistics term which is —1 for bosons and +1 for fermions. 
The energy spectrum has a characteristic Hawking temperature, Th, which is given by 

T« = fti (2.4) 
47rr B H 

and is thus related to the black hole mass and the number of extra dimensions by 

log T H = — log M B h + constant (2.5) 

n + 1 

where the constant is dependent only on n and Mpl- The generator can also model the 
time dependence in which case Th is recalculated after each emission (so as the black hole 
decays it gets hotter). Otherwise the initial X# is used throughout the decay. 
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Due to their high mass, black hole decays are very spectacular events with a large 
visible transverse energy, large multiplicity, and high sphericity with many hard jets and 
leptons. Since most of the black hole decay products result from the evaporation phase, 
as visible Standard Model particles, the ratio of the total hadronic to leptonic activity is 
expected to be roughly 5:1 ||]. A few hard quanta with energy a sizable fraction of the 
Mpl are also expected from the final Planck decay phase ]|]. 

The theoretical work to date has been done in the semi-classical approximation. This 
approximation is only valid if Mbh 3> Mpl 5 -^bh 3* Th and the average multiplicity is 
large, (N) 3> 1. This approximation can only be valid at the LHC if the Planck mass is 
low and even then, there will be problems if the number of dimensions is large (since this 
gives a temperature close to the Planck mass and thus low multiplicity). 



3. Event generation and detector simulation 

CHARYBDIS has been used to generate Monte Carlo event samples. It is interfaced, via the 



Les Houches accord |2^], to HERWIG [23, ^4|] to perform the parton shower evolution of the 



partons produced in the decay and their hadronization. The generated events are then 



passed through the ATLAS fast simulation, ATLFAST [20], in order to give a reasonable 
description of detector resolution and efficiency. 

Unless otherwise mentioned, the CHARYBDIS options were set as follows: 

• Time variation of the black hole temperature was on (TIMVAR=TRUE). 

• Grey-body effects were on (GRYBDY=TRUE) . 

• The black hole was allowed to decay to all Standard Model particles including Higgs 
particles (MSSDEC=3). 

• Kinematic cut-off was turned off (KINCUT=FALSE). 

• The number of particles in the remnant decay was 2 (NB0DY=2). 

This set of options together with the Planck mass set to 1 TeV is called the 'test case' and 
we have used this to illustrate our techniques. Several samples have been generated, so to 
avoid confusion the number of dimensions is always specified. If a mass is given, then the 
generator was forced to produce black holes with a fixed mass, otherwise the range was set 
to 5000-14000 GeV. 

Table |l| summarises the black hole production cross sections at the LHC for n = 2, 4, 
and 6 with M PL = 1 TeV. 2 



4. Model uncertainties 



The theory of black hole production and decay contains many uncertainties and assump- 
tions, particularly at LHC energies. A clear understanding of these is therefore essential 

2 In this analysis, we have used the MRSD-' (DIS) parton distribution function (PDF) set Efj] with 
Q 2 — 1/tbh) where Q 2 is the momentum scale squared at which a PDF is evaluated. 
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Topology 


Total Cross Section (fb) 


5 TeV black hole 


n = 2 


62,000 


n = 4 


37, 000 


n = 6 


34, 000 


8 TeV black hole 


n = 2 


580 


n = 4 


310 


n = 6 


270 


10 TeV black hole 


n = 2 


6.7 


n = 4 


3.4 


n = 6 


2.9 



Table 1: The black hole production cross sections at the LHC for Mpl = 1 TeV as given by 
CHARYBDIS. Note that CHARYBDIS does not include the form factors mentioned in section 0. 



in order for our analyses to be as widely applicable as possible. In this section we review 
these uncertainties. 

4.1 Production cross section 

The process of black hole production in hadron collisions is subject to a number of basic 
uncertainties. The order of magnitude of the parton-level cross section should be given by 



equation 2.1, but the form factor relating the left- and right-hand sides is uncertain and 



would be expected to be n-dependent. Classical numerical simulations |26J suggest values 
in the range 0.5-2, increasing with n. These values are not included in the CHARYBDIS 
generator, but we take them into account when cross section data are used in our analysis 
(in sections and |8|) . 

More fundamentally, the transition from the parton-level to the hadron-level cross 
section is based on the factorization formula 



a(S) = J dxidx 2 f(xi)f(x 2 )o-(s = xix 2 S) (4.1) 



where f(x) is the parton distribution function (PDF) summed over parton flavours. The 
validity of this formula in the trans-Planckian energy region is unclear. Even if factoriza- 
tion remains valid, the extrapolation of the PDFs into this region based on Standard Model 
evolution from present energies is questionable. Also, comparison to Standard Model pro- 
cesses in the trans-Planckian regime would be difficult since perturbative physics would be 
suppressed. 

4.2 The first stages of decay 

CHARYBDIS does not model the initial balding or spin-down phases of the black hole decay. 
The amount of energy emitted from the black hole during these phases is expected to be 
small || so such an omission should not be significant. However, it is probable that the 
energy spectrum will be modified at low energies. 
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4.3 Deposition on the brane 

Estimates vary as to how much energy is expected to be emitted into the bulk via graviton 
emission, but it could be significant. One estimate suggests that the fraction of energy 
emitted into the bulk could be as high as 20% for n = 2-4 rising to nearly 50% for 
n = 7 p7|| . Any energy emitted into the bulk will make an accurate measurement of the 
mass extremely difficult. Although this effect could in principle be observed as a change 
in the expected shape of the cross section as a function of black hole mass, determining 
this would be experimentally challenging. In these studies we have assumed that all the 
energy is deposited on our brane. A modified generator and a more detailed study would 
be necessary to understand the full impact of this assumption. 

4.4 Kinematic limit 

A black hole can only emit a particle with an energy up to half of its mass in order to 
conserve energy-momentum. However, the grey-body distribution used to describe the 
Hawking radiation extends to infinite energy. Although the distribution is only valid for 
very massive (Mbh 3> Mpi,) black holes it is still necessary to deal with the black holes 
as they become lighter. It is expected that given a full theory, the distribution would be 
modified to take this into account, but we have no such theory. 

Figure |l| shows the energy of the primary generator level decay products in the rest 
frame of the black hole. As can be seen, the kinematic limit affects most of the decays. 
This greatly modifies the energy spectrum and also leaves open the question of what to 
do when the generator chooses an unphysical decay, i.e. when it samples from the energy 
spectrum above the kinematic limit. Two options are implemented in CHARYBDIS. In the 
first case (KINECUT=FALSE), if an unphysical decay is chosen, it is thrown away and a new 
one is chosen. This process continues until the black hole has a mass less than the Planck 
mass at which point the decay moves to the final, remnant, stage. In the other option 
(KINECUT=TRUE), when an unphysical decay is chosen, the black hole decay moves straight 
to the final stage. The final stage of the decay is dealt with in the next section. It should be 
noted that for high temperature black holes, where the probability of an unphysical decay 
is large, the different choices implemented in the generator will lead to a large difference 
in the multiplicities and will have a significant impact of the energy distributions. 

4.5 Remnant decay 

At the end of the evaporation phase, a light Planck scale black hole, called a remnant, 
remains which the generator must decay. How this would happen is unknown and will 
only be predicted by a quantum theory of gravity. CHARYBDIS implements this 'remnant' 
decay as an isotropic decay into 2-5 bodies (the number is an option). When the remnant 
decay occurs depends on the option chosen for handling the kinematic limit as described 



in section 4.4. It should be noted therefore, that the uncertainty here can easily affect the 
multiplicity and energy spectra. One example of an affected experimental observable is the 
photon energy spectrum. Figure |2| shows the photon energy distributions (of all photons 
in the event) for 2-body and 4-body remnant decays for two values of n. Even for n = 2 
there is a noticeable effect, but for n = 4, the effect is large. 
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Figure 1: Energy of the generator level decay products in the rest frame of the black hole for a 
5 TeV black hole and 1000 events. The colour scale indicates the number of particles in each bin. 
(a) for n = 2 the kinematic limit (E = Mbh/2, black lines) constricts the energy distribution at 
low masses, (b) for n — 4 the kinematic limit clearly affects the energy distribution at all masses. 
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Figure 2: The photon energy distributions for (a) n = 2 and (b) n = 4. The black and red lines 
are for 2-body and 4-body remnant decays respectively. 



4.6 Time-variation and black hole recoil 

It has been argued || that due to the speed of the decay, the black hole does not have 
enough time to equilibrate between emissions and therefore that the time variation of 
the temperature can be ignored. Therefore, the initial Hawking temperature might be 
measured by fitting Planck's formula for black-body radiation to the energy spectrum of 



the decay products for different bins in the initial black hole mass. Using equation 2.5 the 
number of dimensions can then be extracted. This is the approach taken at a theoretical 
level in ||. 

To illustrate this procedure, we have used the test case with n = 2. Events were 
generated without grey-body factors in 500 GeV mass bins between 5000 and 10000 GeV. 
For each mass bin we have fitted the black-body spectrum to the generator level electron 
energy. Figure ||]a shows the result of this together with the fit using equation |2.5| from 
which we determine n = 1.7 ± 0.3. Figure [3|b shows the result of the same procedure 
and the same test case but with time dependence turned on. In this case we determine 
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n = 3.8 ±1.0 which is well away from the model value. Time dependence is therefore a 
systematic effect with a strong impact on any measurement of n. 




6000 7000 8000 9000 10000 6000 7000 8000 9000 10000 



M BH (GeV) M BH (GeV) 



Figure 3: The plot of Tjj versus Mbh for n = 2 and A/pl = 1 TeV, (a) with a fixed Hawking 
temperature, and (b) with changing (time dependent) Hawking temperature. The statistics used 
correspond to 30 fb" 1 of running at the LHC. 

Another effect that has not been taken into account in previous studies is the recoil 
of the black hole. When a particle is emitted from the black hole, the black hole recoils 
against it. Therefore the next emission is in a boosted frame. Even in the case of a fixed 
temperature decay, the effects of recoil become more significant as the decay progresses 
and the black hole gets lighter. This is exacerbated in the time varying case since the black 
hole also gets hotter as it decays. Any analysis which makes use of the energy spectrum 
should therefore account for this. 

5. Characteristics of the black hole decay 

Black hole decays in the semi-classical limit have high multiplicity. However at LHC 
energies black holes would be on the edge of the semi-classical limit (depending on n) 
which can reduce the multiplicity and make predictions uncertain. This effect can be seen 
in figure |I| which shows that the multiplicity decreases significantly with n. This is due to 
fact that Th is higher for larger n at the same mass. 

A black hole decay is also characterised by a large total transverse energy (figure |||) 
which increases as the black hole mass increases. Even the low multiplicity events tend to 
be rather spherical with high multiplicity events more so. These characteristics are very 
different from standard model and SUSY events which do not have the same access to very 
high energies and tend to produce less spherical events. Therefore, we believe that selecting 
events with high Y^PT, high multiplicity (> 4) and high sphericity will give a pure set of 
black hole events. In addition, it should be noted that the already small Standard Model 
background will be suppressed by the black hole production ||. There are two further 
characteristics which will be interesting to measure and confirm the nature of the events: 
the missing pt {fir) distribution and the charge asymmetry. 
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Figure 4: (a) Average event multiplicities and (b) typical multiplicity distribution with n = 
2, 3, 4, 5, and 6 (black, green, red, cyan, and blue curves, respectively) for 100 fb _1 of integrated 
luminosity. 




Figure 5: The distribution of J^Pt for ro = 2, 3, 4, 5, and 6 (black, green, red, cyan, and blue 
curves, respectively) for 100 fb _1 of integrated luminosity, (a) for 5 TeV and (b) for 8 TeV black 
holes. 



5.1 fx distribution 

Although not all black hole decays contain neutrinos, some will have one or more with 
energies that can be as high as half the black hole mass. The missing energy can be even 
larger than for much of SUSY parameter space. In contrast, most of the Standard Model 
processes tend to have much lower missing transverse momenta. Figure ^ presents the 
distribution of the pT for Standard Model QCD events (with generator level cut px > 
600 GeV), SUSY events (at LHCC SUGRA point 5 @), and 5 TeV black holes with 
n = 2 and 6. 
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5.2 Black hole charge 

Black holes are typically formed from valence quarks, so it is expected that the black 
holes would be charged. The average charge is somewhat energy dependent, but should 
be ~ +2/3. The rest of the charge from the protons is expected to disappear down the 
beam pipes or at very high \tj\. The average black hole charge, (Qbu), can be measured by 
determining the average charge of the charged leptons, (QLept), which should be equal to 
the black hole charge times the probability of emitting a charged lepton. Figure [?] shows 
such a measurement for the test case with n = 2 which gives (Qhept) = 0.1266 ± 0.002 
and thus (Qbu) = 0.654 ± 0.008 using the expected charged lepton emission probability of 
0.1936. 




Figure 6: The distribution of the pr for 
Standard Model QCD events (with generator 
level cut p T > 600 GeV), SUSY events (at 
LHCC SUGRA point 5), and 5 TeV black 
hole with n — 2 and 6. 



Figure 7: The average charge of electrons 
and muons for n = 2 with approximately 
1 ftT 1 of data. 



5.3 Kinematic distributions 

The authors of [29] have studied the hadronic decay of a black hole and found that the 
transverse momentum distribution of charged hadrons depends weakly on the number of 
large extra dimensions. In addition to the event multiplicity and transverse momentum 
distribution, figure ||, we have also looked at the average pr of the events, jets, leptons, and 
the ratio of the difference and sum of the i th and the j th highest p? jet (i ,j = 1, 2, 3, 4) 
and found that these variables also depend only weakly on n. It is therefore not possible 
to get a constraint on n using these distributions. 

5.4 Event shape variables 

In addition to the event multiplicity and spectra, we have studied the following event shape 



variables: the sphericity |30|| , thrust [31], and the Fox- Wolfram moment ratios [g^]. Since 
the sphericity (S) and thrust (T) are sensitive to underlying event and longitudinal motion, 
we have used the corresponding quantities for transverse momenta only. 
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Figure 8: The distribution of (a) the transverse momenta for the highest pr jets, (b) the second 
highest pt jets, (c) the pr of the lcptons, and (d) in the decay of 8 TeV black holes with 
n = 2, 3, 4, 5, and 6 (black, green, red, cyan, and blue curves, respectively) for 100 fb _1 of integrated 
luminosity. 



Defining the Fox- Wolfram moments 

H i = J2 11 p i ( cos %) > 1 = 2 > 3 > - ( 5 - X ) 

i,j vis 

the Fox- Wolfram moment ratios can be expressed as Hi/Hq, where 9ij is the opening angle 
between particles i and j, E v i s is the total visible energy of the event, and P[(x) are the 
Legendre polynomials. 

Figures ^ and [l0| show the distribution of the event shape variables for 5 and 8 TeV 
black holes with n = 2 — 6. The distributions are relatively similar for higher values of 
n, making it hard to distinguish them from one another. Events are relatively spherical 
with a high transverse energy of a few Mpl- For higher dimensions, the events become 
significantly less spherical and are more susceptible to variations in the treatment of the 
remnant decay. 
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(a) (b) 

Figure 9: Event thrust and sphericity for n = 2, 3, 4, 5, and 6 (black, green, red, cyan, and blue 
curves, respectively), (a) The variation of the average thrust and sphericity with mass and n. (b) 
A typical set of distributions of the sphericity for 5 TeV black holes. 
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Figure 10: Fox- Wolfram variables for n = 2, 3, 4, 5, and 6 (black, green, red, cyan, and blue 
curves, respectively), (a) The variation with mass and n. (b) A typical set of distributions for 
5 TeV black holes. 



6. Measurement of the black hole mass 



The black hole 4-momentum is reconstructed simply by summing the 4-momenta of all of 
the particles in the event. We have illustrated this procedure for selected black hole mass 
points 5 and 8 TeV with n = 2 to 6 which were generated with masses ranging 200 GeV 
above and below the selected mass point. Events were selected by requiring at least 4 
jets, all within the acceptance of the ATLAS tracking detector < 2.5). Requiring the 
multiplicity to be greater than 4 ensures that the remnant decay will not be too dominant. 
The transverse momenta of the three highest pt jets was required to be above 500, 400, 
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and 300 GeV respectively. 3 In order to improve the reconstructed mass resolution, events 
were rejected if the missing transverse momentum was greater than 100 GeV. 

The reconstructed Gaussian mass resolution and the overall signal efficiency (the frac- 
tion of accepted events) after the selection cuts for 5 and 8 TeV black hole in n = 2, 4 and 



6 are given in table g with sample plots in figure 11. The mass resolution can be improved 
slightly by raising the threshold of the jet pr, but at the cost of a sharp drop in overall 
signal efficiency. 



Topology 


Mass Resolution (GeV) 


Efficiency (%) 


5 TeV black hole 


n = 2 


202.1 


26.1 


n = 4 


188.4 


30.0 


n = 6 


184.4 


31.9 


8 TeV black hole 


n = 2 


293.9 


13.2 


n = 4 


234.0 


17.8 


n = 6 


226.4 


19.3 



Table 2: The reconstructed Gaussian mass resolution and the overall signal efficiency after the 
selection cuts. 




7. Measurement of the Planck mass 

Some authors || have suggested that since n can be determined from the T^-Mbh rela- 
tionship (equation ^j]), -Mpl can be measured from the normalisation of the temperature. 
For reasons outlined in the next section, we choose not to use this method but instead to 

3 A reconstructed jet was required to have a minimum momentum of 10 GeV within an r\ — 4> cone of 
radius 0.4. 
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follow the suggestion of [^] and determine Mpl from the cross section. In the convention 
used in this paper (and also in []33|]), the cross section is largely independent of n. Figure 12 
shows the parton-level cross section including the corrective form factors calculated in |^6| 
As can be seen, there is very little variation with n. 



■g10 2 



10 




0.5 



1.5 



2.5 
M PL (TeV) 



Figure 12: The partem level cross section for 5 TeV black holes as a function of Mpl for different 
values of n. The form factors calculated in |29l arc included. 



Due to the very high statistics available, the measurement of the parton-level cross 
section will be dominated by the various systematic errors. The main experimental error 
will be the luminosity which should be measured to 5% or better, together with some 
uncertainty in the efficiency This is however likely to be small compared to the theoretical 



uncertainties discussed in section |4.1| . We therefore conservatively estimate that the parton- 
level cross section could be determined to 20% which, for our test case of Mpl = 1 TeV, 
gives an error in Mpl of about 10%. Obviously, the optimal approach is to fit the cross 
section and temperature data simultaneously and this will be demonstrated at the end of 



section 3.2. It is also possible that other processes and observations of new physics at the 
Planck scale may provide independent measurements of Mpl- 



8. Determination of the number of extra dimensions 

Measuring the number of extra dimensions is not a straight-forward task given the uncer- 
tainties outlined in section |2[ One technique that has been suggested Q uses the energy 
spectrum of electrons and photons below Mbh/2. However the authors of || ignore the 
likely effects on the low energy spectrum from the initial parts of the decay (section p~!^ ), 
the effect of the remnant decay (section |4.5| ) and the recoil of the black hole (section [D]) . 
Their analysis is particularly sensitive to these effects because they were attempting to 
use the variation of Tjj with Mbh to measure n. To give some numerical estimates, the 
expected variation, T H (W TeV) - T H (5 TeV) is about 40 GeV for n = 2 and 20 GeV for 
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n = 5 given T#(5 TeV) ~ 200 GeV. This is clearly not a large variation to attempt to 
measure. 

It should also be noted that the differences between different number of dimensions 



become less significant for higher n due to the power law nature of equation 2.5. We have 
investigated a number of variables and techniques whilst considering the effects of the many 
uncertainties. 

8.1 Correlations 

Notwithstanding the comment above, we initially tried to make maximum use of the Mbh~ 
Th relationship. In the case that T# varies with time, there is one Mbh~Th point per 
emission, rather than just one per event. We therefore developed a technique that used all 
the information. The method is: 

1. Reconstruct the black hole from all the particles in the event. 

2. Determine the first(next) particle to be emitted. 

3. Use the measured properties of this particle to determine the temperature of the 
black hole. 

4. Record this mass-temperature point. 

5. Reconstruct the black hole for the next stage using all the particles except for those 
that have been emitted 

6. Repeat steps 2-5 until there are no particles left. 

There are two key parts to this algorithm: determining the order of the emitted particles 
and using the particle properties to determine the temperature. It was hoped that a method 
might be found which when averaged over many events would give the correct Mbb~Th 
relation. 

The method used to determine the order was to assume that the softest particles were 
emitted first. This is because as the black hole decays it gets hotter, and so the average 
energy of the emitted particles should increase. 

The most probable energy and thus pt for the emitted particle is proportional to the 
temperature of the black hole. Since the expected Mbh~Th relation is a power law, the 
Mbb~Pt relation should have the same dependence on n. Thus we have plotted the average 
Pt of the emitted particles for each mass bin. 



Figure 13 shows the result of this method. The shape of the graph is very different 
from that expected, but does show a separation between the numbers of dimensions. This 
is perhaps not surprising since all black hole decays will at some point be affected by the 
kinematic limit problem to which this technique is particularly sensitive. Unfortunately 
this means that extracting the number of dimensions is extremely difficult to do without 
relying purely on Monte Carlo data. Further details on this method are given in [34|. 
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Figure 13: Correlation plots (see text) for 8 TeV black holes (a) with n =2, 3, 4, 5, 6 and (b) with 
M PL =900, 1000, 1100 GeV. 



8.2 Kinematic limit 

In this section we present a new idea that shows a strong variation with the number 
of dimensions and is valid for many different scenarios. It may therefore allow n to be 
determined despite the many uncertainties that affect black hole decays. 

If a particle is emitted with an energy close to the kinematic limit (i.e. E ~ Mbh/2), 
then that particle is probably the first to be emitted. In particular, it is possible to measure 
the fraction of events, p, where the highest energy particle has an energy, E max > E cut 
where E cu t = Mbh/2 — Ed and Ed is a parameter of the analysis that can be chosen, but 
should be small. The probability of the first emission being greater than E cu t can also be 
calculated from integrating the Planck spectrum (equation |2.3| ) thus there is a reasonably 
direct connection between the experimental measurement and theory. 

This method has many advantages: since it deals with the first emission, the probability 
of E max > E cu t will be the same regardless of whether the black hole temperature is time- 
varying or not. Also, since we restrict ourselves to black hole with masses much larger 
than the Planck mass, this measurement will not be affected by Planck scale effects or the 
remnant decay. Any effects which modify the low energy part of the spectrum should also 
have no effect. Finally, the effect of the black hole boost can be taken into account by 
determining E max in the black hole rest frame. 

This technique is however strongly affected by the uncertainty in dealing with the 
kinematic limit (see section since we do not know what the shape of the energy spec- 
trum would be near the kinematic limit. Fortunately this uncertainty can be controlled by 
putting upper and a lower bounds on the theoretical estimate of p. The lower bound is set 
by assuming that none of the unphysical region of the Planck spectrum corresponds to an 
emission with energy greater than E cu t in the actual energy spectrum. The lower bound is 
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therefore 

/•Mbh/2 

Plawer = k / P{E)dE (8.1) 



where A; is a normalising constant, P{E) is the Planck spectrum (equation |2.3j ) and we 
have ignored the effects of the 'grey-body' factors (which are largest at low energy). The 
difference between fermions and bosons is small here. The upper bound uses the opposite 
assumption to the lower bound: that is, that all of the unphysical region in the Planck 
spectrum would correspond to an emission with energy greater than E cu t in the actual 
energy spectrum. The upper bound is therefore 

/"OO 

Pupper = k P(E)dE . (8.2) 

J Ecut 

It should be noted that at very high energies and in the semi-classical limit, p — ► as 
Mbh - ► oo for any fixed Ed- These equations are therefore only valid when p is small. 

There are two competing effects: we would like to set E cu t as low as possible so that 
the upper and lower bound are similar (note that they move apart as E cut — > Mbh/2)- 
However, the lower we set the cut, the greater the chance that the highest emission is 
not the first one. Whilst this is not in itself a problem, it does mean that the differences 
between the time varying and non-time varying cases become more pronounced. In the 
rest of this study we have set Ed = 400 GeV — setting it lower is difficult as experimental 
resolution effects start to dominate. The chance of a soft emission is strongly temperature 
dependent. Points have been omitted if this chance is greater than 50% which affects only 
n = 2, 5 TeV black holes, but would be more significant at lower Mp^. 

We have included an approximate corrective factor for the emission of a first soft 
particle. This was calculated as 

fxE, 

kcor — k 



rxHicut 

\ P(E)dE (8.3) 
J o 



giving p cor = (1 + k cor )p. It is possible for a first emission of any energy to be followed by 
an emission of energy up to the kinematic limit, but the probability of this drops sharply if 
the first emission has an energy greater than 2E cu t. The correction factor, x, was therefore 
set to 2 for the lower bound and 3 for the upper bound. 

A plot of the upper and lower bounds on p for Mpl = 1 TeV is shown in figure 14. As 
can be seen, if Mpl is known, there is good potential for extracting n if it is below 5 and 
could be constrained if n is larger. 

One experimentally tricky aspect of this measurement is that since the black hole 
mass is measured by adding all the particles in the event, the maximum energy in the 
reconstructed black hole rest frame must be less than Mbh/2. This introduces a bias to 
low energies in the measurement of E max and thus reduces the measurement of p. This 
has been corrected for by increasing Ed by 100 GeV which is an estimate for the average 
mis-measurement of E max . It should be noted that not boosting into the black hole rest 
frame would give a very significant over-measurement of p. In addition, to ensure that the 
black hole mass was well measured, events were excluded if they had Pt > 100 GeV, had 
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Figure 14: Upper and lower bounds for p for different n with Mpl = 1 TeV. 



any particles with \r/\ > 2.5 or if the three highest energy particles accounted for more than 
Vo of the total energy in the black hole rest frame. 



Figure 15 shows the results of this analysis for four different models all with n = 4. The 
theoretical upper and lower limits are also shown. For models which have the kinematic cut 
off (see section p~4| ), we would expect them to be consistent with the lower limit; this is the 
case for plots a, b and d in figure |ll| However, if the kinematic limit is on and the remnant 
decay is set to 2-body (plot c), we would expect the data to be consistent with the upper 
limit. This is because if an energy is chosen in the forbidden region (which is the additional 
region included in the upper limit integral), it will cause the decay to be terminated and 
the black hole to split into two. One of the two remnant decay particles must pass the cut. 
As can be seen from the figure, the plots do agree with these expectations. 

Figure 16 shows similar plots to figure 15, but for the test case (see section ||) with 
different values of n. This emphasises that this technique is sensitive to n whilst being 
largely model independent. Note that at high n, the data start to drop below the lower 
limit. This is due to the high temperature here (Tjj ~ 470 GeV for n = 5) which signifi- 
cantly reduces the multiplicity. This result suggests that this analysis has an upper limit 
of validity in the region of Tjj ~ 450-500 GeV. Also, from figure [l6|a, it can be seen that at 
low temperatures, this analysis will not measure the temperature (unless lighter black holes 
are seen, or very high statistics are available), but would instead place an upper limit on 
it. This may be enough to constrain n, or alternatively, measuring the energy distribution 
may be more successful here and these analyses could be combined to give a temperature 
measurement. 

At higher temperatures, this technique becomes strongly sensitive to the temperature 
of the black hole. Indeed the plots above can easily be converted into temperature as a 
function of Mbh- This has been done for the test case with n = 4 in figure |l7] which 
also includes a band equivalent to the upper and lower limits on p. The band includes 
a systematic uncertainty on the measurement of the black hole mass of ±200 GeV. Note 
that as suggested at the beginning of this section, the temperature variation with Mbh 
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Figure 15: Fraction of events passing the cut, p, as a function of Mbh f° r different models all with 
n = 4. Upper and lower bounds for n = 4 are also shown. 



is not determined sufficiently to constrain n without first measuring Mpl. So instead, we 
fix the normalisation of Th at the black hole mass at which it is best measured. In this 
case, we take Tfj = 340 db 30 GeV at Mbh = 7000 GeV. This measurement has been taken 
together with the parton-level cross section with an error of 20% (see section ^) and used 



to determine the model parameters n and Mpl in figure 18. In this case, this gives an 
error on the determination of n of 0.75 and an error on Mpl of 150 GeV. These results are 
indicative of how well this analysis can do. If the cross section error were reduced to 10%, 
the error on Mpl would be 70 GeV and on n, 0.6, showing that, as expected, the cross 
section dominates the determination of Mpl- 

Any improvement in our understanding of how the distribution above the kinematic 
limit should be handled, or how the remnant would decay, would greatly improve this 



analysis by reducing the width of the bands in figures O, \m and 17 



9. Conclusions 

We have discussed the many theoretical uncertainties that can affect black hole decays 
and shown that in at least one case, these can lead to systematic mis-measurements of the 
number of extra dimensions if the analyses previously suggested are used. We have then 
shown the characteristics of black hole decays as they would be measured in the ATLAS 
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Figure 16: Fraction of events passing the cut, p, as a function of Mbh for different values of n for 
the test case. Appropriate upper and lower bounds are shown. 
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Figure 17: Temperature against Mbh for 
n = 4 and 30 fb _1 of integrated luminosity. 
The band shows the systematic uncertainty 
corresponding to the upper and lower bounds 
on p with a systematic on the Mbh measure- 
ment of ±200 GeV. 
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Figure 18: The determination of n and 
AfpL from the measurement of Tjj and an as- 
sumed measurement of the parton-level cross 
section (see text). 



detector. A number of different attempts to determine the model parameters have been 
discussed and a new technique has been introduced. This new technique has been shown 
to control many of the theoretical uncertainties and can be used to measure the black hole 
temperature. We have applied this technique to our test case with four extra dimensions 



and found the temperature to be 340 ± 30 GeV for a black hole mass of 7 TeV. This 
was combined with the parton-level cross section, assumed to be known to 20%, to give 
estimates of the Planck mass and the number of extra dimensions. We conclude that in 
this case the Planck mass can be determined to 15% and the number of extra dimensions 
to ±0.75, with strongly correlated errors. 
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